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ABSTRACT 

Current concepts for Reusable Launch Vehicle design have created renewed interest in 
the use of RP-1 fuels for high pressure and tri-propellant propulsion systems. Such designs 
require the use of analytical methodology which accurately accounts for the effects of real fluid 
properties, combustion of large hydrocarbon fuel molecules, and the possibility of soot 
formation. These effects are inadequately treated in current computational fluid dynamics (CFD) 
codes which are used for propulsion system analyses. 

The objective of this investigation is to provide an accurate analytical description of 
hydrocarbon combustion thermodynamics and kinetics which is sufficiently computationally 
efficient to be practical design tool when used with CFD codes such as the FDNS code. 

A rigorous description of real fluid properties for RP-1 and its combustion products will 
be derived from the literature and from experiments conducted in this investigation. Upon the 
establishment of such a description, the fluid description will be simplified by using the 
minimum of empiricism necessary to maintain accurate combustion analyses and including such 
empirical models into an appropriate CFD code. An additional benefit from this approach is that 
the real fluid properties analysis simplifies the introduction of the effects of droplet sprays into 
the combustion model. 

Typical species compositions of RP-1 have been identified, surrogate fuels have been 
established for analyses, and combustion and sooting reaction kinetics models have been 
developed. Methods for predicting the necessary real fluid properties have been developed and 
essential experiments have have been designed. Verification studies are in progress, and 
preliminary results from these studies will be presented. The approach has been determined to 
be feasible, and upon its completion the required methodology for accurate performance and heat 
transfer CFD analyses for high pressure, tripropellant propulsion systems will be available. 
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RP-1 COMBUSTION CHEMISTRY 


• RP-l combustion model w/soot formation 

• Verify model w/data from literature 

• Verify model w/new test data 

• Describe real fluid thermodynamic properties 

• Add real fluid HC properties to CICM 

• Real fluid single element model 

Real fluid HC, H 2 , and 0 2 
Real fluid tri-propellant 
Additional turbulence models 

• Account for radiation 

• Assemble elements to make injector model 
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RP-1 


RP-1 is a straight run kerosene fraction which is subjected to acid washing and sulfur 
dioxide extraction. The mean molecular weight is about 174, and the H/C ratio is about 1.9. 
This implies that RP-1 is a multicomponent hydrocarbon fuel without a specified species 
distribution. 
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Distillation Curve for RP-1 
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| HYDROCARBONS ISOLATED FROM ONE REFERENCE PETROLEUM 

Molecular 

Formula 

Name 

Type 

Normal 
Boiling Point 

C 10 H 14 

1.3-Diethylbenzene 

Benzene 

181.10 

c 10 h 14 

l-Methyl-3-propylbenzene 

Benzene 

181.80 

Ck)H 14 

n-Bulylbenzene 

Benzene 

183.27 

C 10 H 18 

Bicycloparaffin 

Bicycloparaffin 

183.4 

CiqHis 

Bicycloparaffin 

Bicycloparaffin 

183.7 

C I0 H U 

1 -Methyl-4-propylbenzene 

Benzene 

183.30 

c l 0 H 14 

l, 2-Diethylbenzene 

Benzene 

183.42 

C,oH 14 

1 ,3-Dimethyl-5-ethylbenzene 

Benzene 

183.58 

C 10 H 14 

1, 4-Diethylbenzene 

Benzene 

183.78 

CioH 14 

1, Methyl-2-propylbenzene 

Benzene 

184.80 

C l0 H 14 

1, 4-Dimethyl-2 -ethylbenzene 

Benzene 

186.83 

C 10 H 18 

trans-Decahydro-naphthalene 

Bicycloparaffin 

187.25 

C 10 H 14 

1, 3-Dimethyl-4-ethylbenzene 

Benzene 

188.20 

CiqH | 4 

1, 2-Dimethyl-4-ethylbenzene 

Benzene 

189.48 

C 10 H 16 

Tricyclo (3. 3. 1.1)- decane 

Tricycloparaffin 

190 

c 10 h 14 

1, 3-Dimethyl-2 -ethylbenzene 

Benzene 

190.01 

C 10 H 12 

1-Methylindan 

Aromatic- 

cycloparaffin 

190.6 

c 10 h 12 

2-Methyl indan 

Aromatic- 

cycloparaffin 

191.4 

C 10 H l4 

1, 2-Dimethyl-3 ethylbenzene 

Benzene 

193.91 

C 10 H l8 

cis-Decahydro-naphthalene 

Bicycloparaffin 

195.69 

C u H24 

n-Undecane 

Normal paraffin 

195.89 

C 10 H 14 

1,2,4,5-Tetramethyl-benzene 

Benzene 

196.80 

Cio^14 

1 ,2 , 3 , 5-Tetr amethyl -benzene 

Benzene 

198.00 

C11H20 

Bicycloparaffin 

Bicycloparaffin 

202.5 

C„H 16 

1 -Methyl-3-n-butylbenzene 

Benzene 

204.1 


1 ,2,3,4-Tetramethyl-benzene 

Benzene 

205.04 


4-Methyl indan Aromatic- 

cycloparaffin 


CiqHu 
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HYDROCARBONS ISOLATED FROM ONE REFERENCE PETROLEUM (Continued) 


Molecular 

Formula 

Name 

Type 

Normal 
Boiling Point 

c„h 16 

l,3-Dimethyl-4-n propylbenzene 

Benzene 

206.6 

C 10 H 12 

1,2,3,4-Tetra-hydronaphthalene 

Aromatic-cycloparaffin 

207.57 

C„H 16 

1 ,2-Dimethyl-4-n-propylbenzene 

Benzene 

208.5 

C„H l6 

Trimethylethylbenzene 

Benzene 

212.3 

^12^26 

n-Dodecane 

Normal paraffin 

216.28 

CjoHg 

Naphthalene 

Dinuclear aromatic 

217.96 

C„H 14 

2-Methyl-(l,2,3,4-tetra- 

hydronaphthalene) 

Aromatic-cycloparaffin 

220.7 

c„h 14 

6-Methyl-(l,2,3,4-tetra- 

hydronaphthalene) 

Aromatic-cycloparaffin 

229.03 

c 14^50 

2,6, 10-Trimethylundecane 

Branched paraffin 

231 

c„H I4 

5-Methyl-(l,2,3,4-tetra- 

hydronaphthalene) 

Aromatic-cycloparaffin 

234.35 


n-Tridecane 

Normal paraffin 

235.43 

CjiHjo 

2-Methylnaphthalene 

Dinuclear aromatic 

241.05 

C„H 10 

1-Methylnaphthalene 

Dinuclear aromatic 

244.64 

C 15 H 32 

2,6,10-Trimethyldodecane 

Branched paraffin 

249 


n-Tetradecane 

Normal paraffin 

253.52 

c 12 h 10 

Biphenyl 

Dinuclear aromatic 

255.0 

Cj 3 H| 2 

2-Methylbiphenyl 

Dinuclear aromatic 

255.3 

C 12 H 12 

2-Ethylnaphthalene 

Dinuclear aromatic 

257.9 

C 12 H 12 

1-Ethylnaphthalene 

Dinuclear aromatic 

258.7 

^12^12 

2 ,6-D imethylnaphthal ene 

Dinuclear aromatic 

262 

C 12 H 12 

2,7-Dimethylnaphthalene 

Dinuclear aromatic 

263 

Ci 2 H 12 

1 ,7-Dimethylnaphthalene 

Dinuclear aromatic 

263 

C, 2 H 12 

1 ,6-Dimethylnaphthalene 

Dinuclear aromatic 

263 

C 12 H 12 

1 ,3-Dimethylnaphthalene 

Dinuclear aromatic 

265 

C 12 H l2 

1 ,5-Dimethylnaphthalene 

Dinuclear aromatic 

265 

c 14 h 14 

2,5-Dimethylbiphenyl 

Dinuclear aromatic 

267 

C,5H 16 

Trimethylbiphenyl 

Dinuclear aromatic 

267 

C 16 Hi 8 

T etramethy lb iphenyl 

Dinuclear aromatic 

267 
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RP-1 SURROGATE FUELS 




^Mol % [ 

MB 

C 13 H 12 

methylbiphenyl 

17.4 

255 

c 12^24 

n-heptylcyclopentane 

45.4 

224 

C| 2 H 2 g 

n-tridecane 

37.2 

235 



100.0 



M wt = 173.9 
H/C = 1.922 


Note: Critical pressure for RP-1 is 340 psia, 
and critical temperature is 679°K. 
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RP-1 COMBUSTION PROPERTIES 


HEAT OF COMBUSTION 

HOC = -18640 (Btu/lbm) = -10.346 (kcal/gm) 


HEAT OF VAPORIZATION 

HOV = 106 (Btu/lbm) 


HEAT OF FORMATION (HOF) 

To determine a HOF for mixtures, an effective molecular formula must be specified and 
used to evaluate the HOF. Frequently, an arbitrary molecular weight of 100 gms is assumed 
for a basis for thermodynamic calculations. 


H/C 

HOF(liccai/ibOgm) 

Source 


1.8624 

-36.01 

Lockheed 

-10.287 

1.90 

-41.6 

SAIC 

-10.305 

2.0 

-44.36 

Aerojet ’78 

-10.403 

1.9063 

-42.0 

TMX - 1783 

-10.305 

1.9423 

-33.068(v) 

SP-273 

-10.441(v) 

1.9423 

-38.946(1) 

SP-273 

-10.382(1) 

1.922 

-39.915 

Surrogate Fuel 

-10.346 
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Compound 


SPECTROSCOPIC ANALYSIS 
OF H-l GAS GENERATOR EXHAUST GASES 


Sample #1 Sample #3 Sample #4 


Propyne 


Propene 


Diacetylene 


1,3 Butadiene 


2 Butene 


- Butene 


1.5 Hexadyne 


3 Methyl Pentene -1 


Cyclopentene 


1,2,3 trimethylcyclopantane 


Ethyl Benzene 


Carbon 


RP-1 


26.15* 


38.51 


17.58 


37.93 



17.14 
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Generalized Combustion Kinetics Model 
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Initial Polycyclic Aromatic Hydrocarbon (PAH) Formation 
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Generated needed C p data in CEC format 

S°298 corrected for symmetry and optical isomers using Benson’s 
data. 



Aerosol Dynamics Based on Frenklach and Harris’s “Method II 
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Cp’s for Diphenyl Compounds 
Dotted lines are synthesized data 



y/do 
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Toluene and Iso-Octane Quasiglobal Kinetics Model 
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SOOT PARAMETERS 
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[soot] = mass of soot/volume of gas (g/cm 3 ). 



Elementary Reactions 
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RP-1 COMBUSTION VERIFICATION 
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SOOT CONCENTRATION, mg/1 iter on dry basis 



EQUIVALENCE RATIO 


Comparison Of Quasiglobal Model Prediction Of Soot Emissions 
With Experimental Data For Toluene And Toluene/Iso-Octane 
Blends. Computation: Well-Stirred Reactor; Experiment: 
Jet-Stired Combustor. 
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NET SOOT EMISSION («9/l ) 



JSC Soot Emissions For Toluene/Methylnaphthalene Blends; 
T = 1900 K; r = 3 ms. 
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Soot Concentration (ng/1 9 300K, I <tn) 



Dependence Of Soot Concentration On Equivalence Ratio; 
Ratio; Atomization Air = 15 g/min/nozzle; Indicat ed 
Temperature = 1900 ±35K. 
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852 


Photograph of the exit flow at an exit pressure of 2 atm. Region 1 is the undisturbed cone 
region 2 is the Prandtl-Meyer expansion zone, and region 3 is bhe mixing zone. The 
spectrometer line of sight i6 indicated by a dotted line. 
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Wavelength (/Am) 
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Wavelength (/im) 
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Wavelength 





| v. ;**>•. vff Soot Density from GDC Experiments 

1HBBI 

HBUi 

lilllll 

wmMi 

« x (cm 2 /gin) 

0. 1400E-02 

1680.0000 

1.8500 

0. 1087E-06 

0. 1288E+05 

0.7000E-02 

1415.0000 

1.6000 

0.5200E-06 

0.1346E+05 

0.3200E-01 

1290.0000 

1.5000 

0.2343E-05 

0. 1366E+05 

0.9000E-01 

1045.0000 

1.3300 

0.653 IE-05 

0. 1378E+05 

0. 1000E-01 

1760.0000 

1.6000 

0.7886E-06 

0. 1268E+05 

0.5700E-01 

1470.0000 

1.4500 

0.4268E-05 

0.1335E+05 

0. 1850E+00 

1230.0000 

1.2500 

0. 1348E-04 

0.1373E+05 

0. 1050E-01 

2600.0000 

2.1000 

0.9944E-06 

0.1056E+05 

0.2550E-01 

2370.0000 

1.8500 

0.2295E-05 

0.1111E+05 

0.1250E+00 

1900.0000 



EB9 

0.3900E+00 

1700.0000 

1.4000 

0.3040E-04 

0.1283E+05 


Rayleigh Theory for Small Particles 


K y 


36k nk 

K {n 2 -k 2 +l) 2 + 4 n 2 k 2 



K x (cm'‘) = /c x p c = linear absorption coefficient 
* x (cm 2 /gm) = mass absorption coefficient 
P c (gm/cc of volume) = soot density 
p 0 (gm/cc of soot) = bulk soot density 
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